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A previous publication [Europhysics Letters 78, 10009 (2007)] suggested to coherently generate 
mesoscopic superpositions of a two-component Bose-Einstein condensate in a double well under 
perfectly symmetric conditions. However, already tiny asymmetries can destroy the entanglement 
properties of the ground state. Nevertheless, even under more realistic conditions, the scheme is 
demonstrated numerically to generate mesoscopic superpositions. 
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I. INTRODUCTION 



Generation of mesoscopic entanglement in a Bose- 
Einstein condensate (BEC) [J, i, i, H i, H 0, is 
still a challenge of fundamental research. A possible sys- 
tem for which entanglement could be achieved is a BEC 
in a double- well potential [H, Q. While attractive con- 
densates can become unstable [loj . for repulsive single 
species condensates the ground state is not the desired 
superposition (see, e.g., Ref 0]). Nevertheless, by phase 
engineering fidelities (i.e. the probability to be in an ideal 
superposition) of about 60% have been obtained numer- 
ically [|. 

In Ref. [ll[ we suggested to use a binary condensate to 
produce Bell-like mesoscopic superpositions with fideli- 
ties above 95% by using a method of tunneling control 
which was very recently realized experimentally in an op- 
tical lattice [12 ] . However, to achieve the entanglement 
generation, the ground state of a two-component con- 
densate of exactly the same number of particles of each 
kind was followed adiabatically in a perfectly symmetric 
double well potential [l3j] ■ Given the fact that even very 
small asymmetries can destroy entanglement properties 
of eigen-functions [3] , this rises the question if the gen- 
eration of mesoscopic entanglement proposed in Ref. [ill ] 
would also be possible under more realistic conditions, 
i.e., slightly different particle numbers and a small tilt. 

The focus of this manuscript is to investigate the 
scheme of entanglement generation including asymme- 
tries in both the potential and the particle numbers. In 
Sec.[TT]the model used to describe a two-component BEC 
in a double well is introduced. While Sec. [TTT] explains 
how the mesoscopic superpositions can be characterized, 
Sec. IIVI outlines the noticeable differences between the 
ideal case considered in Ref. UM and the more realistic 
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FIG. 1: (color online) Sketch of one of the target states for 
entanglement generation in a double well: a Bell-like super- 
position of two wave-functions for a binary Bose-Einstein con- 
densate. 



(slightly asymmetric) situation considered here for which 
the entanglement generation no longer is entirely adia- 
batic. 



II. MODEL, CONTROLLING THE 
ENTANGLEMENT GENERATION 

To describe a binary Bose-Einstein condensate in a 
double well, we use the Hamiltonian in two-mode ap- 
proximation (Ref. [HI cf. Ref. [l6|): 



H = — — [a ia ' 2 + a[a 2 
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(1) 



where (b^) creates, respectively annihilates, a bo- 
son of type A (B) in well j, Q is the tunneling frequency, 
ka (^b) is the interaction parameter for atoms of type 
A (B) and kab models the interaction between distin- 
guishable particles. Such Hamiltonians have also been 
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used to describe entanglement generation beyond BEC- 
physics [l7j . The (unwanted) small tilt which might not 
be avoidable experimentally is represented by e. 

The time-dependent potential difference oc f(t) can 
be used to perform tunneling control via resonant tun- 
neling @, El] (cf. Ref. 0). References @, [IH generate 
multi-particle entanglement by periodic driving, 



f(t) = sin(wt) 



(2) 



to effectively renormalize the tunneling frequency Q by 
the Jo-Bessel function in the high frequency limit to 



fieff = SlJo(2,u/ w ) 



(3) 



This mechanism for tunneling control (cf. p^|). which 
on the single particle level has been investigated theo- 
retically already for some 16 years [2(1 Hi], [22], HH, has 
very recently been observed experimentally for a Bosc- 
Einstein condensate in an optical lattice [l2j . 

Another possibility to achieve the entanglement gen- 
eration of Ref. [ll| is to decrease the tunneling fre- 
quency by increasing the distance d between the two 
wells. Approximating the wave-functions of particles of 
mass m in the two wells by the ground states of har- 
monic oscillators of oscillation frequency ui and oscillator 
length £ osc = y / h/(2muj), one obtains jj, HH 



exp 



i 2 

v osc 



for the tunneling frequency and 

Ha 



(4) 



(5) 



where a is the s-wave scattering length, for the interaction 
parameters of Eq. |T]) as a function of Q,. 

A possible realization of such a model could be two 
hyperfme-states of, e.g, 87 Rb. Feshbach resonances for 
interaction between different hyper-fine states have been 
observed experimentally [24| . thus providing a way to 
slightly increase the interaction between distinguishable 
particles as compared to identical particles. 

To coherently generate entanglement, Ref. 



Hi sug- 



gests to decrease the (effective) tunneling frequency for 
such a situation, thus adiabatically following the ground 
state towards highly entangled mesoscopic states. As 
decoherence-times tend to rapidly decrease with increas- 
ing particle numbers, not too large particle numbers 
should be used to investigate entanglement generation. 
Experimentally, Bose-Einstein condensates of the order 
of 100 particles have been realized [2^|; in Ref. [ll[ we 
used a binary BEC with exactly Na = Nb = 50 parti- 
cles of each kind. Here, deviations Na = 50 — AN/2, 
N B = 50 + AN/2 up to \AN\/(N A + N B ) = 10% will be 
considered. 
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FIG. 2: Fidelity ((7} as a function of dimensionless time r = fit 
if the tunneling frequency is reduced linearly from fi(r = 0) = 
1 to Q(t = 500) = and then held at this value. Solid line: 
p+, dashed line p_. The double well is slightly asymmetric 
(s = 0.001) with iV A = 45, N B = 55, k a = k b = 0.04; 

KAB = 0.05. 

Inset: The same situation for an untilted double well. Even 
without the tilt, the ratio kab/ka has to be larger than in 
Ref. [Tl| ] in order to ensure that the Bell-state still is the 
ground state for a symmetric potential. 



III. IDENTIFYING MESOSCOPIC 
SUPERPOSITIONS 

For the symmetric conditions of Ref. [ll| the fidelity 
obtained by projecting the wave-function \ip) on the Bell 
state 
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(\0,N a }a\N b ,0)b± |JV a ,0)a|0,JV, 



B B 



(6) 



provides an excellent measure of entanglement genera- 
tion, 



p± = \m±)\ 2 



(7) 



if P± > 0-5- In Fig. [2] the tunneling frequency in an 
asymmetric double well containing a binary BEC with 
Na 7^ N-g is reduced linearly during half the time and 
during the other half £1 held at zero: 



n = 



i - 
o 



< T < T, 



ramp 



= m 



(8) 



For a slightly tilted double well, the ground state can- 
not be one of the two superpositions of Eq. ©: In the 
limit fl — > the two wave-functions used in the super- 
positions ^ do have different energies. Thus, even if a 
perfect superposition was reached, the phase would os- 
cillate similar to Fig. [5J The fact that the fidelities stay 
below 1 indicates that more Fock states then just those 
included in Eq. Increasing the tilt by, say, a factor 
of ten would further increase the oscillation frequency as 
the energy difference between both wells increases while 
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increasing it by a much larger factor would destroy the 
generation of mesoscopic superpositions. Periodic driv- 
ing introduces an additional frequency, thus leading to 
even larger oscillations frequencies in the fidelities. 

For the purpose of identifying mesoscopic superpo- 
sitions independent of these oscillations, we define the 
modified fidelity: 

1 2 
Pfid = ^ (VPi + VP2) > ( 9 ) 

where 

Pl = {(y\0,N A ) A \N B ,0) B f (10) 

and 

P2 = «*|iVA,0)A|0,7V B )B) 2 • (11) 

As for the usual fidelity, this measure is 1 for a perfect 
superposition and already p^d w 0.5 might signify that 
there is no superposition at all. The modified fidelity (|9|) 
can be derived by projecting the wavefunction on the 
superposition 

IW> - ^(|0,JV A )|JV B ,0) +e*\N A ,0)\0,N B )) (12) 

and maximizing the resulting projection 

^ = l(%)f (13) 

as a function of ip. 

IV. ENTANGLEMENT GENERATION ONLY 
PARTIALLY ADIABATIC 

In order to understand why generation of mesoscopic 
superpositions is still possible, Fig. [3] shows the mod- 
ified fidelity © for different ramping times r ramp (see 
Eq. ©). If the ramping is done on too short time- 
scales, the oscillations characteristic for too fast, i.e. non- 
adiabatic, ramping occur and high fidelities cannot be 
reached. If, on the other hand, the ramping is done too 
slowly, the system does not reach a mesoscopic superpo- 
sitions either as it adiabatically follows the ground state 
for too long and thus becomes rather asymmetric. For 
moderate ramping, the system first adiabatically follows 
the ground state thus leading to a wave-function with a 
bimodal probability distribution. Then, the symmetry 
of the wave-function is conserved (instead of tunneling 
through energetically higher states with half the parti- 
cles of at least one kind in both of the wells) to finally 
reach a mesoscopic superposition. 

If the tunneling control is done via slowly variing the 
amplitude fx of periodic driving as in Ref. [11| (cf. Eqs. (0) 
and |3])), i.e., £1 = const., cu ^> and 

2.4048w t 

V - 5 ■ < r < r ramp , (14) 



T3 




FIG. 3: Modified fidelity j9) as a function of time t. For 
ramping time r ra mp = 1000 (solid line), pad reaches 95.61% 
by first adiabatically following the ground state and then non- 
adiabatically conserving the symmetry of the wave-function. 
The other parameters are chosen as in Fig [5] A ramping time 
of T ramp = 50 (dotted line) is much too short as can be seen 
in the oscillations indicating non-adiabaticity. On the other 
hand T ramp = 50000 (dashed line) is too slow; it adiabatically 
follows the ground state up to values where it becomes rather 
asymmetric. This leads to a final state with a modified fi- 
delity well below 90%. 

Inset: For a ramping time of r ramp = 1000, the final modified 
fidelity is shown as a function of the particle number differ- 
ence A with Na = 50 — A/2 and iV B = 50 + A/2. Up to 
A/50 = 20%, the value used for most of the figures, a trans- 
fer to highly entangled states is achieved. For higher A, the 
ramping-time and/or the ratio of the interaction parameters 
would have to be adjusted. 



where 2.4048 ... is the first zero of the Jo-Bessel function 
(see Eq. (|3])), the generation of mesoscopic entanglement 
also occurs. For the parameters in Fig. [5] with a ramping 
time of r r amp = 1000, Fig. [4] shows the emergence of an 
entangled wave-function from an initial state for which 
the probability distribution is not bimodal. 



V. CONCLUSIONS 

While already two independent Bose-Einstein conden- 
sates provide interesting effects (see, e.g. Ref. [26j]), two- 
component BECs provide an even wider variety of ef- 
fects [53, [28[. Here it has been shown that for a bi- 
nary condensate in a double-well potential, the genera- 
tion of mesoscopic entanglement suggested in Ref. [ll[ 
works also under less ideal conditions, i.e., if the parti- 
cle numbers are different and if the double well is not 
perfectly symmetric. 

For a given experimental situation, the linear ramping 
used here in the spirit of coherent control theory can be 
optimized both toward short ramping times and large 
fidelities, e.g., by methods of optimal control theory [29j |. 
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FIG. 4: (color online) The projection of the wave-function on 
Fock-states \i, Na — i)\j, Nb — j) for a binary condensate for 
which the tunneling control is performed via periodic driving 
as suggested in Ref. [ll| (see Eqs. (fl4|> and (0). While the 
initial state at r = (top) is no bimodal distribution at all, 
a mesoscopic superposition is generated already at r = 400 
(bottom). The parameters where chosen as in Fig [3] with a 



ramping time of T ri 



1000. The driving frequency to 
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the final fidelity is again well above 90%. 



We would like to thank Y. Castin for his support. N. T. 
thanks the Institut Henri Poincare-Centre Emilc Borcl 
for hospitality and support. Funding by the European 
Union (C. W. through contract MEIF-CT-2006-038407 
and N. T. through contract MEST-CT-2005-019755) is 
gratefully acknowledged. 



[1] J. I. Cirac, M. Lewenstein, K. M0lmer, and P. Zoller, 
Phys. Rev. A 57, 1208 (1998). 

[2] A. S0rensen, L. M. Duan, J. I. Cirac, and P. Zoller, Na- 
ture 409, 63 (2001). 

[3] M. Holthaus and S. Stenholm, Europhys. J. B 20, 451 
(2001). 

[4] V. I. Yukalov, E. P. Yukalova, and V. S. Bagnato, Phys. 

Rev. A 66, 043602 (2002). 
[5] A. Micheli, D. Jaksch, J. I. Cirac, and P. Zoller, Phys. 

Rev. A 67, 013607 (2003). 
[6] K. W. Mahmud, H. Perry, and W. P. Reinhardt, J. Phys. 

B 36, L265 (2003). 
[7] C. Weiss and T. Jinasundera, Phys. Rev. A 72, 053626 

(2005). 

[8] C. E. Creffield, Quantum control and entanglement us- 
ing periodic driving fields (2007), arXiv: 0704 . 1792 vl 



[cond-mat . other] . 
[9] R. Gati and M. K. Oberthaler, J. Phys. B 40, R61 (2007). 

[10] R. J. Dodd, M. Edwards, C. J. Williams, C. W. Clark, 
M. J. Holland, P. A. Ruprecht, and K. Burnett, Phys. 
Rev. A 54, 661 (1996). 

[11] N. Teichmann and C. Weiss, EPL 78, 10009 (2007). 

[12] H. Lignier, C. Sias, D. Ciampini, Y. Singh, A. Zen- 
esini, O. Morsch, and E. Arimondo, Dynamical control 
of matter-wave tunneling in periodic potentials (2007), 
larXiv : 070"T 0403vl [cond-mat . other] . 

[13] As the coherent control was performed via time-periodic 
potential differences, strictly speaking the corresponding 
Floquet-state was followed adiabatically, see also [30( ] . 

[14] D. R. Dounas-Frazer and L. D. Carr, Tunneling reso- 
nances and entanglement dynamics of cold bosons in the 
double well (2006), [arXiviquant^p h/OeiOieelrl. 



■5 



[15] G. J. Milburn, J. Corney, E. M. Wright, and D. F. Walls, 

Phys. Rev. A 55, 4318 (1997). 
[16] H. J. Lipkin, N. Meshkov, and A. J. Click, Nucl. Phys. 

62, 188 (1965). 

[17] S. Dusuel and J. Vidal, Phys. Rev. B 71, 224420 (2005). 

[18] C. Weiss, Las. Phys. Lett. 3, 212 (2006). 

[19] S. Haroche, C. Cohen- Tannoudji, C. Audoin, and J. P. 

Schermann, Phys. Rev. Lett. 24, 861 (1970). 
[20] F. Grossmann, T. Dittrich, P. Jung, and P. Hanggi, Phys. 

Rev. Lett. 67, 516 (1991). 
[21] M. Holthaus, Phys. Rev. Lett. 69, 1596 (1992). 
[22] M. Grifoni and P. Hanggi, Phys. Rep. 304, 229 (1998). 
[23] C. Weiss, Phys. Rev. B 73, 054301 (2006). 
[24] A. Marte, T. Volz, J. Schuster, S. Diirr, G. Rempe, 

E. G. M. van Kempen, and B. J. Verhaar, Phys. Rev. 



Lett. 89, 283202 (2002). 
[25] C.-S. Chuu, F. Schreck, T. P. Meyrath, J. L. Hanssen, 
G. N. Price, and M. G. Raizen, Phys. Rev. Lett. 95, 
260403 (2005). 

[26] Y. Castin and J. Dalibard, Phys. Rev. A 55, 4330 (1997). 
[27] C. J. Myatt, E. A. Burt, R. W. Christ, E. A. Cornell, 

and C. E. Wieman, Phys. Rev. Lett. 78, 586 (1997). 
[28] R. J. Ballagh, K. Burnett, and T. F. Scott, Phys. Rev. 

Lett. 78, 1607 (1997). 
[29] A. Assion, T. Baumert, M. Bergt, T. Brixner, B. Kiefer, 

V. Seyfried, M. Strehle, and G. Gerber, Science 282, 919 

(1998). 

[30] H. P. Breuer and M. Holthaus, Z. Phys. D 11, 1 (1989). 



